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This study aimed to explore the neural development status of the visual system of children (around 8 years 
old) using contrast sensitivity. We achieved this by eliminating the influence of higher order aberrations 
(HOAs) with adaptive optics correction. We measured HOAs, modulation transfer functions (MTFs) and 
contrast sensitivity functions (CSFs) of six children and five adults with both corrected and uncorrected 
HOAs. We found that when HOAs were corrected, children and adults both showed improvements in MTF 
and CSF. However, the CSF of children was still lower than the adult level, indicating the difference in 
contrast sensitivity between groups cannot be explained by differences in optical factors. Further study 
showed that the difference between the groups also could not be explained by differences in non-visual 
factors. With these results we concluded that the neural systems underlying vision in children of around 8 
years old are still immature in contrast sensitivity. 

The contrast sensitivity function (CSF) shows the relationship between the scale of a stimulus and its 
minimum detectable contrast. Measuring the CSF is one of the most fundamental assessments of visual 
function. The shape of the CSF is jointly determined by the optical and neural components of the visual 
system^"^. Previous studies show worse contrast sensitivity in children, and suggest this mainly stems from an 
immature neural system^"^. Children greater than 4 years old typically no longer have the refractive and accom- 
modative errors that are common during infancy, therefore these optical factors will not affect the performance of 
these children^. 

This old story needs to be updated however due to developments in our knowledge of the optical system, 
especially higher-order aberrations (HOAs)^ of the eye. HOAs are aberrations beside refractive errors (defocus 
and astigmatism), which cannot be measured or corrected by traditional methods^. Several studies have shown 
that HOAs vary among individuals ages^^'^^, pupil sizes^^'^^, ethnicities refractive errors^^'^^ and even 
time^°. The hallmark effect of HOAs is that it reduces the resolution of images of the retina^\ Contrast sensitivity 
has been demonstrated to increase when HOAs are corrected with an adaptive optics (AO) system in normal 
adults^'^^'^^ and in the elderly^^. Larger values of HOAs have been found in children than in adults in several 
studies^^'^^'^^, however their effects on the contrast sensitivity of children are still not fully understood. 

Since HOAs vary among individuals and affect the measurement of contrast sensitivity, a precise comparison 
that eliminates the influence of HOAs is necessary for exploring the status of children's visual neural system in 
contrast sensitivity. In this study, we measured contrast sensitivities at different spatial frequencies in children of 
around 8 years old and in young adults under both monochromatic HOAs-corrected and HOAs-uncorrected 
conditions. Based on the measurements, we addressed one question: whether the neural component of the visual 
system of children of around 8 years old has reached an adult level in terms of contrast sensitivity? The reason for 
choosing children of around 8 years old is that previous studies have not reached an agreement on the devel- 
opmental status of 8 -year- old children in contrast sensitivity^"^'^^"^\ 

We found that after HOAs were corrected both children and adults showed improvements in contrast sens- 
itivity, while adults still exhibited a better performance than children. The worse performance of children may 
result from an immature visual neural system or simply non-visual factors (e.g., head movements, eye movements 
and inattention etc.). These non-visual factors would increase measured threshold^^'^^ and within -subject vari- 
ability^'^^'^^. To reduce their influence, we conducted a 10-day training in children. After training, the within- 
subject variability of the CSF test in children decreased to the same level as adults, which suggested that these 
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non-visual differences may have little influence on the CSF compar- 
ison between the two groups in the post-training measurement. In 
this condition, adults still exhibited a better performance than chil- 
dren. Our results indicate a substantial difference in spatial contrast 
sensitivity between children of around 8 years old and adults, which 
provides a new insight into the development of spatial vision in 
human. 

Results 

In general, subjects (all have normal vision without correction) had 
tests of contrast sensitivity function (CSF) conducted monocularly 
with a randomly chosen eye in the dark room on a real-time closed- 
loop adaptive optics (AO) visual stimulator system^^. Two CSFs were 
tested on different days: one of them was measured with higher- order 
aberrations (HOAs) of the eye being corrected in real-time during 
the test ("with AO correction" in our figures) and the other one was 
measured with HOAs being uncorrected during the test ("without 
AO correction"). Optical aberrations were recorded during the CSF 
tests. Contrast sensitivity was defined as the reciprocal of contrast 
threshold for detecting sine-wave gratings at 79.3% correct. 

Correction of optical aberrations. Figure 1 illustrates the average 
root mean square (RMS) value of Zernike coefficients up to 7th order 
with and without AO correction. Naturally, the children had larger 
HOAs (t(9) = 2.436, p = 0.038, two-tailed Independent Samples t- 
test) and third-order aberrations (t(9) = 3.157, p = 0.012) than the 
adults, but similar total aberrations {p > 0.05) and low-order 
aberrations (astigmatism and defocus) as adults (p > 0.05). The 
AO system worked well in correcting optical aberrations, which 
reduced the average RMS value of wavefront aberrations from 
0.202 ± 0.016 |im (mean ± SE) to 0.085 ± 0.004 |im for the child 
group and from 0.187 ± 0.019 jim to 0.071 ± 0.002 jim for the adult 
group. Both reductions were significant (for children, t(5) = 6.634, p 
= 0.001; and for adults, t(4) = 6.209, p = 0.003). 

We also compared coma-like and spherical aberrations in the two 
age groups both with and without AO correction, as previous reports 
have shown that these aberrations could have a large influence on the 
contrast sensitivity function^^"^^. We found these aberrations were 
not significantly different between groups {p > 0.05) either with or 
without AO correction. 

We calculated the horizontal (as the orientation of the sine-wave 
grating we used was vertical) modulation transfer function (MTF)^° 
to further quantify the improvement of the optical system by AO 
correction. The MTF indicates the ratio of image contrast on the 
retina to object contrast as a function of spatial frequency, which is 
a straightforward measure of the image quality on the retina^\ The 
average MTFs with and without AO correction in the two age groups 
are presented in Figure 2a (to compare with the contrast sensitivity 



functions, we focused on modulation transfer factors of the seven 
spatial frequencies that were tested in the contrast sensitivity mea- 
surement). We found comparable MTFs in the two groups (F(l,9) = 
0.064, p = 0.806) when HOAs were not corrected. AO correction 
significantly improved the MTFs in both the child (F( 1,5) = 39.872,p 
= 0.001) and adult (F(l,4) = 26.798, p = 0.007) groups. After AO 
correction, the MTFs in the two groups still overlapped and had a 
small but significant difference (F(l,9) = 8.290, p = 0.018) due to 
small between- subject variance in MTF. 

Contrast sensitivity. The average contrast sensitivity functions 
(CSFs) with and without AO correction in the two age groups are 
plotted in Figure 2b. When HOAs were not corrected, the contrast 
sensitivity of the child group was significantly lower than that of the 
adult group (F(l,9) = 24.108, p = 0.001). The interaction between 
spatial frequency and age group was significant (F(6,54) = 18.641, 
< 0.001): larger differences were shown at lower spatial frequencies. 
When HOAs were corrected, the CSFs were improved in both 
children (F(l,5) = 7.498, p = 0.041) and adults (but not signi- 
ficantly, F(l,4) = 3.235, p = 0.146). However, the CSFs of children 
were still worse than those of adults (F(l,9) = 15.761, p = 0.003), 
even worse than the performance of adults that were measured 
without AO correction (F(l,9) = 12.943, p = 0.006). Similar to the 
situation when HOAs were not corrected, the interaction between 
spatial frequency and age group was also significant when HOAs 
were corrected (F(6,54) = 8.921, p < 0.001), with larger diffe- 
rences found at lower spatial frequencies. 

Performance of the children after 10 days of training. We have 
shown that our child subjects still exhibit lower contrast sensitivity 
than adults even when HOAs were corrected. The worse perfor- 
mance of children could result from immature neural circuits 
underlying spatial contrast sensitivity, but could also be due to the 
influence of non-visual factors (such as head movements, eye 
movements, inattention etc.). Since the non-visual factors have 
been shown to increase both measured threshold^^'^^ and within- 
subject variability^'^^'^^ we examined the extent of these non- visual 
effects by comparing the within-subject variability of the CSF test 
between the child and the adult groups. We found that the within- 
subject variability of the CSF test with AO correction (i.e., standard 
error of the reversal points of contrast in each staircase) was 
significantly larger in children than in adults (Figure 3, F(l,9) = 
7.296, p = 0.024). To rule out the possible influence of non-visual 
factors on our measured thresholds we conducted 10 days of training 
for the children. A high spatial frequency, at which the contrast 
detection threshold was 0.4 in pre-training CSF test with AO 
correction, was selected as training frequency. The same task as 
CSF measurement (i.e., contrast detection task) was used for the 
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Figure 1 | Average root mean square (RMS) value of Zernike coefficients with and without adaptive optics (AO) correction. Average RMS value of 
wavefront aberrations was provided up to 7th order. In all conditions, un-filled bars represent results that were measured without AO correction; filled 
bars represent results that were measured with AO correction; blue bars represent results of children {n = 6) and red bars represent results of adults {n = 
5). Error bars represent standard errors of RMS value across subjects. represents p < 0.05 by a two-tailed Independent Samples t-test. 
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Figure 2 | Effect of AO correction on Modulation transfer functions (MTFs) and Contrast sensitivity functions (CSFs) in the child and the adult 
groups, (a) Modulation transfer factor with and without AO correction as a function of spatial frequency for the child and the adult groups. Modulation 
transfer factor of the eye denotes the contrast ratio of one particular spatial frequency that remains after the image comes through the eye; the x-axis shows 
spatial frequency (cycles/°). The black solid Une represents the diffraction limit. Error bars indicate standard errors of the modulation transfer factors 
across subjects, (b) Contrast sensitivity with and without AO correction as a function of spatial frequency for the child and adult groups. Contrast 
sensitivity was defined as the reciprocal of contrast threshold for detecting sine-wave gratings at 79.3% correct. Error bars show standard errors of the 
contrast sensitivity across subjects. 



training. Training successfully reduced the within- subject variability 
of the post-training CSF measurement in children such that it was no 
longer different from that in adults (Figure 3, F(l,9) < 0.001, p = 
0.993). However, even though the contrast sensitivity of children 
with AO correction tended to improve after training (Figure 4a, 
F(l,5) = 6.201, p = 0.055), it was still worse than that of adults 
(Figure 4b, F(l,9) = 9.466, p = 0.013). The interaction between 
spatial frequency and age group was also significant (F(6,54) = 
5.332, p < 0.001), and larger differences were found at lower 
spatial frequencies (Figure 4b). 

To further investigate the differences in contrast sensitivity 
between the trained children and the adults, we also conducted a 
curve-fitting approach to compare the CSFs between the two groups. 
We fitted the CSFs with difference of Gaussian (DOG) functions and 
used an F-test for nested models^^ to compare a full model (where 
different DOG functions are needed to fit the two CSFs) with a 
reduced model (where the same DOG function is used to fit the 
two CSFs). The full model generated a quite good fit to our data 
(R2 _ 99 5%)^ which was significantly better than the reduced model 
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Figure 3 | Average within- subject variability of the CSF test of the 
untrained and trained children and untrained adults. The within- subject 
variability in the CSF test was defined as the standard error of the reversal 
points (contrast in dB) in each staircase. The average within- subject 
variability of the CSF test was averaged across subjects and spatial 
frequencies. The blue bar shows results from the untrained children, the 
purple bar shows results from the trained children, and the red bar shows 
results from the adults. represents p < 0.05 by a two-tailed Independent 
Samples t-test. Error bars represent standard errors of the variability of the 
CSF test across subjects. 



(F(3,4) = 17.521, p = 0.009). These results suggested that two dif- 
ferent DOG functions were needed to fit CSFs of the trained children 
and the adults. The fitted results based on the full model are shown in 
Figure 4b as solid lines. Obviously, the CSF of children is lower than 
that of adults. 

Discussion 

We found that AO correction successfully improved optical quality 
by decreasing the wavefront aberrations and increasing the modu- 
lation transfer factor in both the child and the adult groups. AO 
correction also improved the CSF in both groups. Under the 
HO As corrected condition, the contrast sensitivity of the child group 
was still significantly lower than that of the adult group. After intens- 
ive training, the within-subject variability of the CSF test with AO 
correction in children was reduced to a comparable level as in adults. 
Therefore in this condition the non-visual differences may have little 
influence on any differences between the two groups. However, the 
contrast sensitivity of the children in the trained condition was still 
worse than that of the adults. These results indicate an immaturity in 
the neural circuits underlying spatial contrast sensitivity in children. 

One technical concern is the influence of optical factors such as 
fluctuation of accommodation^^ and low order aberrations (defocus 
and astigmatism) ^'^^'^^ on the measurement of the MTF and CSF. For 
the fluctuation of accommodation^^ we corrected subjects' aberra- 
tions in real time during the whole CSF test. In order to reduce the 
low order aberrations' effects, all subjects we chose had normal vision 
without spectacles. The RMS value of the wavefront aberrations also 
showed that there was no difference between children and adults on 
the second order aberrations. On the other hand, we used a fixed 
defocus correction to correct chromatic aberration for both children 
and adults. It might or might not be the best position for all subjects, 
but from their subjective reports and their C-letter-test results we 
were confident that it was still in their clear position. Nevertheless, we 
should point that some potential uncorrected defocus might still exist 
and thus affect subjects' performance. However, such uncorrected 
defocus, if existed, was likely to have no notable effect on the com- 
parison of CSFs between children and adults in the current study: 
Firstly, it has been shown that defocus has larger effect on visual 
performance with AO-correction than that without AO-correc- 
tion^^'^^'^^. Thus, if the uncorrected defocus affected children more, 
we would expect less benefit from AO correction in children, which 
was not the case in our study. Secondly, several studies have demon- 
strated that defocus mainly damage contrast sensitivity at higher 
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Figure 4 | Contrast sensitivity functions (CSFs) with AO correction in trained and untrained children and in untrained aduhs. (a) CSFs with AO 
correction in the trained ('# ') and untrained ('A') children. Error bars represent the standard error in contrast sensitivity across subjects, (b) CSFs with 
AO correction in the trained children ('# ') and adults ('■'). Error bars represent the standard error in contrast sensitivity across subjects. The solid lines 
are curve-fits with difference of Gaussian (DOG) fijnctions. 



spatial frequencies, but not at low spatial frequencies^'^^'^^. While 
when compared the CSF of children with adults, larger differences 
were found at lower spatial frequencies but not at higher frequencies. 

Another technical concern is the influence of non-visual factors 
such as head movements, eye movements and inattention on the 
comparison between adults and children. Head movements and 
eye movements were carefully controlled in our study. During the 
test, subjects were instructed to keep stable with their heads on a chin 
rest to minimize the head movements. Subjects were also asked to 
maintain their fixation at the center of the screen to reduce the eye 
movements. On the other hand, stimulus was shown immediately 
after fixation and its presentation time was too short (i.e., 117 ms) for 
subjects to saccade^". Effects of inattention were reduced by asking 
subjects to take a break for 2-3 minutes after fmishing the test in one 
block (typically lasted for about 5 minutes). Though we tried our best 
to minimize the infiuence of these non-visual factors, they still may 
have affected our comparison between children and adults, since the 
variability of our CSF test was significantly larger in children than in 
adults. To further diminish any such effects, we conducted a 10-day 
training protocol with a contrast detection task at a high spatial 
frequency. Training successfully reduced the variability of the CSF 
test of children to the same level as adults. In this condition therefore 
these non-visual factors may have a comparable infiuence on chil- 
dren and adults, however the contrast sensitivity of children was still 
worse than that of adults. It should be noted that former studies had 
shown that training with the contrast detection task would bring 
about perceptual learning in adult and child amblyopes^^"^^ and in 
normal adults under AO condition^^. It is quite possible that our 
training in children also induced some perceptual learning effect, 
which in turn underestimated the difference of contrast sensitivities 
between children and adults. Therefore, we may confidently con- 
clude that the contrast sensitivity of children is lower than that of 
adults. 

Our data showed that HOAs were higher in children than adults, 
which is consistent with several previous reports^^'^^'^^ Meanwhile, 
total aberrations and MTFs in the child and adult groups were similar 
before AO correction. The similarity of MTF in children and adults 
in the current study indicates similar optical quality in these two 
groups. This is consistent with a previous study which showed that 
the average MTF of children was slightly better than that of adults, 
but there was considerable overlap between MTFs for the two age 
groups^^. 

The aim of our study is to check whether the contrast sensitivity of 
the visual neural system in children of around 8 years old reaches the 
adult level. We showed that the CSF of children with AO correction 



was lower than that of adults, which can't be attributed to optical and 
non-visual factors. One notable phenomenon is that the peak con- 
trast sensitivity of our subjects was not high even in adults (—100 in 
adults). Such peak contrast sensitivity is lower than that in some 
reports^'^'^^'^^'^^, but also consistent with the magnitude in several 
other studies^^'^^'^^'^\ The variance of peak contrast sensitivity in 
different studies might come from different paradigms that were 
used in testing. For example, our study used a quite short exposure 
time of stimuli (117 ms) during the measurement, while most other 
studies presented stimuli for several seconds^'^^ or with unlimited 
time^^'^^ in testing the contrast sensitivity of children and adults. 
Nevertheless, we believe that the magnitude of the peak contrast 
sensitivity couldn't limit the ability to detect the differences between 
children and adults in contrast sensitivity, as there were several stud- 
ies that reported higher peak contrast sensitivity but still found less 
contrast sensitivity in children of around 8 years old than that of 
adults^^'^\ 

Our data, therefore support the existence of immature neural 
circuits underlying spatial contrast sensitivity in children of around 
8 years old^^. This idea is consistent with several neurophysiological 
studies^^"^^. For instance, experiment on macaque monkeys sug- 
gested that the main factors limiting spatial contrast sensitivity in 
development are in the nervous system, not in the eye or photore- 
ceptors^^. Further study showed that the asynchronous changes also 
exist between the LGN cells and the behavioral performance^^, indi- 
cating a later maturation in the visual cortex than in the LGN. The 
immature contrast sensitivity of our child subjects indicates the exist- 
ence of a subsequent process of neural maturation in the visual 
system after the age of 8 years old. 

Methods 

Subjects. Six children (Mean age ( ± SD): 8.5 ± 0.55 years old; 2 males) and five adults 
(25 ± 2.00 years old; 3 males) with normal vision (Snellen acuity equal to or better 
than 20/20 in each eye without wearing any optical correction) participated in this 
experiment. All subjects were naive to the purpose of the experiment. A written 
informed consent was obtained from each adult subject and each child's guardians/ 
parents before the experiment. This study was approved by the Institutional Review 
Board of University of Science and Technology of China and the Institute of Optics 
and Electronics in the Chinese Academy of Sciences. The methods were carried out in 
accordance with the approved guidelines. 

Apparatus. All the measures were conducted on a real-time closed-loop adaptive 
optics (AO) visual stimulator system^*^. The display had a pixel resolution of 800 X 
600 and a frame refresh rate of 60 Hz. A special circuit was used to produce 14 bits of 
gray levels^^. 

Subjects viewed a 1.5° test field monocularly through a 4-mm artificial pupil in a 
dark room. Their pupils were naturally dilated after several minutes' adaptation in the 
dark room before the measurements were made. Their heads were fixed with a chin 
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rest to rule out any influence of head movements during the test. They were instructed 
to focus on the center of the display when performing the task. Retinal illuminance, 
defined as the product of the luminance in the pupil plane (18.9 cd/m^) with the pupil 
area, was 237 Troland. 

The minidisplay was fixed on a stepping motor for adjusting the potential focus 
shifts. Two experienced adult operators performed the adjustment of the stimulus 
position for a perfect focus for a Landolt C letter with all ocular aberrations corrected 
for several times. The position of the stepping motor was set to the average of these 
repetitions. In pilot study, we confirmed that all our subjects could clearly see the C 
letter (better than 20/20) at this position. Those who could not see clearly at this 
position were excluded in our experiment. 

Aberrations calculation. Aberrations were calculated up to 35 Zernike coefficients 
(7th order). Root mean square (RMS) value of wavefront aberrations was calculated 
as the square root of the sum of the squares of the corresponding coefficients from 
each Zernike expansion. The modulation transfer function (MTF, describing the 
attenuation of contrast as a function of spatial frequency of the optical system) was 
calculated from the wavefront aberrations based on the standard Fourier optics 
theories*^". Since vertical sine-wave gratings were used in our experiment, only the 
horizontal MTF was calculated. Contrast sensitivity was measured at seven spatial 
frequencies (0.9, 1.5, 3, 6, 12, 24, and 36 cycles/"). 

Experiment design. Before the measurement, all subjects practiced for one hour in 
order to understand the experiment and familiarize themselves with the 
manipulations. 

RMS values, MTFs and CSFs were then measured under both the HOAs-corrected 
(i.e. with AO correction) and HOAs-uncorrected (i.e. without AO correction) con- 
ditions. AO correction operated in real-time along the whole CSF measurement 
under the HOAs-corrected condition. CSF measurements in the two conditions were 
finished separately in two different days with counterbalanced orders in different 
individuals. 

To further diminish the influence of non-visual factors, a 10-day training protocol 
was conducted after the measurements. A high spatial frequency, at which the con- 
trast detection threshold was 0.4 in pre-training CSF test with AO correction, was 
selected as the training frequency. The same contrast detection task as was used in the 
CSF measurement was used for the training. After training, CSFs were re-measured 
under the HOAs-corrected condition using the same paradigm as the pre-training 
test. 

Procedure. A two -interval forced-choice contrast detection task was used for both 
training and CSF measurements. The presentation sequence in each trial was as 
follows: a 267-ms fixation signaled by a brief sound, a 117-ms stimulus interval 
(such a short period of presentation was used to minimize the influence of eye 
movements^"), a 500-ms blank, a 267-ms fixation signaled by a brief sound, a 117- 
ms second stimulus interval, a 500-ms blank, and blank until response. A vertical 
sine-wave grating was randomly presented in one of the two stimulus intervals, 
while the other stimulus interval was blank. Subjects were instructed to answer 
which stimulus interval contained the signal by pressing a key. The next trial 
started immediately after the response. The contrast of the sine-wave grating was 
varied by a three-down one-up staircase procedure with a proportional step size of 
10% to track the 79.3% correct value. For CSF measurement, each session 
consisted of 8 blocks with each containing 77 trials. Different spatial frequencies 
were randomly intermixed in each block. For training, one session consisted of 7 
blocks with each containing 90 trials, and a fixed spatial frequency (the frequency 
at which the contrast threshold was 0.4 in the pre-training CSF test with AO 
correction). Subjects were instructed to take a break for 2-3 minutes after 
finishing one block. The CSF measurement and training in each day lasted about 
50 minutes in total. 

Data analysis. The contrast detection threshold at each spatial frequency was 
calculated by averaging the contrast of the reversal points in the corresponding 
staircase. 

The within-subject variability of the CSF test was defined as the standard error of 
the reversal points (contrast in dB) in each staircase. Contrast in dB was calculated as 
20 X logio (contrast). The average within- subject variability of the CSF test was 
averaged within-subject variability across subjects and spatial frequencies. 

The differences of RMS values between the two age groups were compared using 
two-tailed Independent Samples t-test. The differences of MTF, CSF and variability 
between the two groups were compared using repeated- measures ANOVA with one 
within- subjects factor (spatial frequency) and one between- subjects factor (age 
group). The differences of RMS values with and without AO correction in each group 
were compared using two-tailed Paired Samples t-test. The differences of MTF or CSF 
with and without AO correction in each group were compared using repeated- 
measures ANOVA with two within- subjects factors (spatial frequency and AO 
correction). 

The method of curve-fitting for CSFs is similar to a previously documented method 
in Huang et aV\ In particular, the average CSFs of children and adults were fitted 
with both a full model and a reduced model. For the full model, the CSF of child group 
was fitted with a difference of Gaussian (DOG) functions (G1-G2). The CSF of adult 
group was fitted with (G1-G2 + G3). The model assumed that the difference between 
the two age groups was a Gaussian function: G3. For the reduced model, the CSFs of 
children and adults were both fitted with a same function: G1-G2. This model 



assumed that there was no difference between the two groups. We used an F-test for 
nested models"^^ to compare these two models. 



{^full ^reduced] 
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Where dfi = kjuu — Keducedy and df2 = N — Keduced — hk represents the number of 
parameters in each model and N is the number of data points. When p < 0.05, we 
chose the full model, if not, we chose the reduced model. 
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